99 FILES FILO,FILO
100 PRINT "&&&L100&&&STARTING"
105 SCRATCH #1
110 READ D$,O$
111 WRITE #1,D$
112 WRITE #1,O$
120 READ S1,S2,S3,N,F1,F2
121 WRITE #1,S1
122 WRITE #1,S2
123 WRITE #1,S3
124 WRITE #1,N
125 WRITE #1,F1
126 WRITE #1,F2
127 PRINT "---L127---";S1;S2;S3;N;F1;F2
130 LET S4=S2+S3
140 LET S5=S4+1
150 READ M,V1
151 WRITE #1,M
152 WRITE #1,V1
160 FOR I=1 TO S4
170 READ S(I)
171 WRITE #1,S(I)
180 NEXT I
190 FOR I=1 TO S4
200 READ A$(I)
201 WRITE #1,A$(I)
210 NEXT I
220 FOR I=1 TON 
230 READ L(I),R(I)
231 WRITE #1,L(I)
232 WRITE #1,R(I)
240NEXT I
250 READ T1,T2,T3,T4
251 WRITE #1,T1
252 WRITE #1,T2
253 WRITE #1,T3
254 WRITE #1,T4
260 LET E2=0
262 LET E3=0
270 LET E2=FNT(T2)-FNT(T1)
280 IF S1=1 THEN 300
290 LET E3=FNT(T4)-FNT(T3)
300 LET S6=S4+2
305 PRINT S1;S2;S3;S4;S5;S6
310 FOR I=1 TO N
320 READ P(I)
321 WRITE #1,P(I)
330 NEXT I
400 GOSUB 1000
500 GOSUB 2000
550 IF L$="PILOT1" THEN 400
600 GOTO 99999
1000 REM THIS IS ORIGINALLY PILOT1
1010 DIM A$(18),S(18),L(100),P(100),R(100),Y(100)
1020 DIM C(20,20),Z(20,20)   
1030 DIM Q(20),X(20)
1033 DIM D(100),N(100)
1035 PRINT "-----L1035 PILOT1----"
1038 RESTORE #1
1040 READ #1,D$,O$
1050 READ #1,S1,S2,S3,N,F1,F2
1060 LET S4 = S2 + S3
1070 LET S5=S4+1
1080 READ #1,M,V1
1090 FOR I = 1 TO S4 
1100 READ #1,S(I)
1110 NEXT I
1120 FOR I = 1 TO S4 
1130 READ #1,A$(I)
1140 NEXT I
1150 FOR I = 1 TO N
1160 READ #1,L(I),R(I)
1170 NEXT I
1180 READ #1,T1,T2,T3,T4
1190 LET E2=E3=0   
1200 LET E2 = FNT(T2)-FNT(T1)
1210 IF S1=1 THEN 1230
1220 LET E3 = FNT(T4)-FNT(T3)
1230 LET S6=S4+2   
1235 PRINT S1;S2;S3;S4;S5;S6
1240 FOR I = 1 TO N
1250 READ #1,P(I)
1260 NEXT I
1270 MAT Z = ZER(S6,S6)      
1280 MAT C = ZER(S6,S6)      
1290 FOR I=1 TO S6 
1300 LET X(I)=Q(I)=O 
1310 NEXT I
1320 FOR I = 1 TO N
1330 LET Y(I) = F1*(P(I)-F2) 
1340 IF S1=1 THEN 1430
1350 IF L(I)>S2 THEN 1380
1360 LET Y(I) = Y(I)-E2      
1370 GO TO 1390
1380 LET Y(I) = Y(I)-E3      
1390 IF R(I)>S2 THEN 1420
1400 LET Y(I) = Y(I)+E2      
1410 GO TO 1430
1420 LET Y(I) = Y(I)+E3      
1430 LET L = L(I)  
1440 LET R = R(I)  
1450 LET Q(L) = Q(L)+Y(I)    
1460 LET Q(R) = Q(R)-Y(I)    
1470 LET Z(L,L) = Z(L,L)+1   
1480 LET Z(R,R) = Z(R,R)+1   
1490 LET Z(L,R) = Z(L,R)-1   
1500 LET Z(R,L) = Z(R,L)-1   
1510 LET Z(S5,L) = Z(S5,L)+1 
1520 LETZ(L,S5) =Z(S5,L)     
1530 LET Z(S5,R) = Z(S5,R)-1 
1540 LET Z(R,S5) = Z(S5,R)   
1550 LET Z(S5,S5) = Z(S5,S5)+1 
1560 LET Q(S5) = Q(S5)+Y(I)  
1570 NEXT I
1580 FOR I = 1 TO S4 
1590  LET Z(S6,I) = S(I)     
1600 LET Z(I,S6) = S(I)      
1610 NEXT I
1620 MAT C = INV(Z)
1630 FOR I = 1 TO S6 
1640 FOR J = 1 TO S6 
1650 LET X(I) = X(I) +C(I,J)*Q(J)      
1660 NEXT J
1670 NEXT I
1680 SCRATCH#1     
1690 WRITE #1,D$
1691 WRITE #1,O$
1700 WRITE #1,S1
1701 WRITE #1,S2
1702 WRITE #1,S3
1703 WRITE #1,S4
1704 WRITE #1,S5
1705 WRITE #1,S6
1710 WRITE #1,T1
1711 WRITE #1,T2
1712 WRITE #1,T3
1713 WRITE #1,T4
1720 WRITE #1,N
1721 WRITE #1,M
1722 WRITE #1,V1
1723 WRITE #1,F1
1724 WRITE #1,F2
1725 WRITE #1,E2
1726 WRITE #1,E3
1730 FOR I=1 TO S4 
1740 WRITE #1,A$(I)
1741 WRITE #1,S(I)
1750 NEXT I
1760 FOR I=1 TO N  
1770 WRITE #1,L(I)
1771 WRITE #1,P(I)
1772 WRITE #1,R(I)
1773 WRITE #1,Y(I)
1780 NEXT I
1790 FOR I=1 TO S6 
1800 WRITE #1,X(I)
1810 NEXT I
1820 DEF FNT(T)=(T-20)*(-40.6+(T-20)*(-.95+.01*(T-20)))    
1830 PRINT 
1840 PRINT 
1850 PRINT "*****THISSHOULD BE A CHAIN*****"
1855 RETURN
2000 PRINT
2010 PRINT
2020 PRINT
2030 PRINT TAB(25);"S A N D I A"
2040 PRINT
2050 PRINT"          LEAST SQUARE ANALYSIS OF STANDARD CELL DATA"     
2060 PRINT
2070 PRINT 
2080 PRINT 
2090 PRINT 
2120 RESTORE #1
2130 READ #1,D$
2131 READ #1 O$
2140 READ #1, S1, S2, S3, S4, S5, S6   
2150 READ #1, T1, T2, T3, T4 
2160 READ  #1, N,M,V1,F1,F2,E2,E3      
2170 FOR I=1TO S4  
2180 READ #1, A$(I), S(I)    
2190 NEXT I
2200 FOR I=1 TO N  
2210 READ #1, L(I),P(I),R(I),Y(I)      
2220 NEXT I
2230 FOR I=1TO S6  
2240 READ #1, X(I) 
2250 NEXT I
2260 LET V2 = 0    
2270 FOR I = 1 TO N
2280 LET L = L(I)  
2290 LET R = R(I)  
2300 LET D(I) = Y(I)-(X(L)-X(R)+X(S5)) 
2310 LET V2 = D(I)*D(I)+V2   
2320 NEXT I
2330 LET V3=N-S4   
2340 LET W3=0      
2350 LET V4=SQR(V2/V3)       
2360 LET W2 = 0    
2370 PRINT "DATE-----";D$;TAB(30);"OBSERVER-----";O$       
2380 PRINT" DESIGN----";S2;"X";S3;",",N;"OBSERVATIONS"     
2390 PRINT 
2400 PRINT " STD DEV ="; V4,"ACC STD DEV ="; V1     
2410 PRINT " DEG OF FREEDOM ="; V3      
2420 PRINT 
2430 LET W1 = 0    
2440 IF W2 = 0 THEN 2460
2450 GO SUB 3640
2460 FOR I = 1 TO N 
2470 LET L = L(I) 
2480 LET R = R(I) 
2490 IF D(I)=0 THEN 2520
2500 LET X1=ABS(D(I))       
2510 LET D(I)=(D(I)/X1)*(INT(X1*100+.5))/100    
2520 IF W2 = 1 THEN 2560
2530 IF ABS(D(I))<3*V1 THEN 2580
2540 LET W1 = W1 + 1
2550 IF W1 = 1 THEN 2600
2560 PRINT I;TAB(6);A$(L);TAB(16);A$(R);TAB(27);
2561 PRINT P(I);TAB(35);Y(I);
2570 PRINT TAB(49); D(I)     
2580 NEXT I       
2590 GO TO 2630
2600 PRINT"       THE FOLLOWING DEVIATIONS EXCEED 3 SIGMA"
2610 GO SUB 3640
2620 GO TO 2560
2630 IF W2 = 1 THEN 2650
2640 IF W1 = 0 THEN 2670
2650 GO SUB 3670
2660 GO TO 2680
2670 PRINT" NO DEVIATIONS EXCEED 3 SIGMA."      
2680 IF W2 = 1 THEN 2740
2690 PRINT" FOR A COMPLETE LISTING OF DEV. ENTER 1, OTHERWISE ENTER 0"
2700 INPUT W2     
2710 PRINT
2720 IF W2 = 1 THEN 2430
2730 IF W1 = 0 THEN 2800
2740 PRINT" ENTER NO. OF OBS. TO BE DELETED"    
2750 INPUT N5     
2760 IF N5 <> O THEN 3260
2770 PRINT
2780 PRINT
2790 PRINT
2800 PRINT
2810 PRINT
2820 IF M <>0 THEN 2850
2830 PRINT" MEAN OF GROUP WAS NOT GIVEN."       
2840 GO TO 2860
2850 PRINT "REF GRP MEAN REDUCED BY 1.01, IN MICROVOLTS IS";M       
2860 PRINT
2870 PRINT" LEFT-RIGHT COMPONENT =";X(S5)       
2880 PRINT" TEMP GRP1 ="; T2,"TEMP CORR GRP1 ="; E2 
2890 IF S1 = 1 THEN 2910
2900 PRINT" TEMP GRP2 =";T4,"TEMP CORR GRP2 =";E3 
2910 PRINT
2920 PRINT "GRP1 CORR TO";T1;"DEG"    
2930 IF S1=1 THEN 2960
2940 PRINT
2950 PRINT "GRP2 CORR TO";T3;"DEG"    
2960 PRINT
2970 PRINT" ESTIMATES OF DIFF FROM MEAN AND TEMP CORRECTED EMFS"    
2980 PRINT
2990 PRINT "CELL NO"," ","EST VALUE","EMF 101+" 
3000 GO SUB 3670
3010 LET P5=0     
3020 FOR I = 1 TO S4
3030 PRINT A$(I), 
3040 IF S(I) = 1 THEN 3090
3050 PRINT" UNK ",
3060 LET P5=P5+1  
3070 LET W3=W3+X(I) 
3080 GO TO 3100
3090 PRINT" STD ",
3100 PRINT X(I),  
3110 IF M = 0 THEN 3140
3120 PRINT INT((X(I)+M)*100+.5)/100   
3130 GO TO 3150
3140 PRINT" "     
3150 NEXT I       
3160 IF M = 0 THEN 3220
3170 LET W3=W3/P5 
3180 LET W3 = INT((W3+M)*10+.5)       
3190 GO SUB 3670
3200 PRINT
3210 PRINT" MEAN OF UNK GRP REDUCED BY 1.01 IS";W3
3220 FOR I = 1 TO 5 
3230 PRINT
3240 NEXT I       
3250 STOP 
3260 PRINT N5;"LINE(S) ARE TO BE DELETED"       
3270 PRINT
3280 LET H=0      
3290 FOR K=1 TO N5
3300 PRINT "ENTER ONE OBSERVATION NUMBER TO BE DELETED"   
3310  INPUT N(K)  
3320 PRINT
3330 NEXT K       
3340 FOR I=1 TO (N5-1)      
3350 FOR J=1 TO (N5-I)      
3360 IF N(J)<N(J+1) THEN 3400
3370 LET T5=N(J)  
3380 LET N(J)=N(J+1)
3390 LET N(J+1)=T5
3400 NEXT J       
3410 NEXT I       
3420 PRINT "THE FOLLOWING OBSERVATIONS WILL BE DELETED"   
3430 PRINT
3440 FOR K=1 TO N5
3450 LET N6=N(K)  
3460 LET L = L(N6)
3470 LET R = R(N6)
3480 PRINT N6;TAB(6);A$(L);TAB(16);A$(R);TAB(27);P(N6)    
3490 NEXT K       
3500 PRINT
3510 FOR K=1 TO N5
3520 LET N6=N(K)  
3530 LET N6=N6-H  
3540 LET N = N-1  
3550 FOR I = N6 TO N
3560 LET P(I) = P(I+1)      
3570 LET L(I) = L(I+1)      
3580 LET R(I) = R(I+1)      
3590 NEXT I       
3600 LET H=H+1    
3610 NEXT K       
3620 PRINT
3630GOTO 3690
3640 PRINT" OBS";TAB(7);"LEFT";TAB(17);"RIGHT";TAB(27);"OBS.";       
3650 PRINT TAB(38);"OBS";TAB(49);"DEV."
3660 PRINT" NO.";TAB(7);"CELL";TAB(17);"CELL";TAB(38);"T-CORR"      
3670 PRINT"----------------------------------------------------------"
3675 LET L$="THATS ALL"
3680 RETURN       
3690 SCRATCH #1   
3700 WRITE #1,D$
3701 WRITE #1,O$
3710 WRITE #1,S1
3711 WRITE #1,S2
3712 WRITE #1,S3
3713 WRITE #1,N
3714 WRITE #1,F1
3715 WRITE #1,F2
3720 WRITE #1,M
3721 WRITE #1,V1
3730 FOR I=1 TO S4
3740 WRITE #1, S(I) 
3750 NEXT I       
3760 FOR I=1 TO S4
3770 WRITE #1, A$(I)
3780 NEXT I       
3790 FOR I= 1 TO N
3791 WRITE #1,L(I)
3792 WRITE #1,R(I)
3793 NEXT I
3800 WRITE #1,T1
3801 WRITE #1,T2
3802 WRITE #1,T3
3803 WRITE #1,T4
3830 FOR I = 1 TO N
3840 WRITE #1, P(I) 
3850 NEXT I       
3860 PRINT
3865 LET L$="PILOT1"
3866 PRINT "GOING BACK TO PILOT1"
3870 RETURN
9000 DATA JAN-70
9001 DATA M8/GU-RUN1
9030 DATA 2,6,6
9045 DATA 24,1,0
9060 DATA 8240.70,1
9075 DATA 1,1,1,1,1,1,0,0,0,0,0,0
9090 DATA G1-1,G1-2,G1-3,G1-4,G1-5,G1-6
9105 DATA M8-1,M8-2,M8-3,M8-4,M8-5,M8-6
9120 DATA 7,2,7,4,11,4,11,6,9,6,9,2
9135 DATA 8,1,8,3,12,3,12,5,10,5,10,1
9150 DATA 1,11,1,9,5,9,5,7,3,7,3,11
9165 DATA 2,12,4,12,4,8,6,8,6,10,2,10
9180 DATA 28.,28.001,28.,28.001
9195 DATA 2.7,5.3,6.,2.1,12.8,4.0
9210 DATA 15.8,7.4,5.7,12.6,12.6,14.1
9225 DATA -14.2,-14.9,-13.4,-14.2,-5.0,-5.7
9240 DATA -3.5,-5.9,-7.6,-13.8,-12.1,-3.3
99999 END
                                                                                                                                                                                                                                                                                                                                  